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We study the radiative heat transfer and the Casimir-Lifshitz force occurring between two bodies 
in a system out of thermal equilibrium. We consider bodies of arbitrary shape and dielectric prop- 
erties, held at two different temperatures, and immersed in a environmental radiation at a third 
different temperature. We derive explicit closed-form analytic expressions for the correlations of the 
electromagnetic field, and for the heat transfer and Casimir-Lifshitz force, in terms of the bodies 
scattering matrices. We then consider some particular cases which we investigate in detail: the 
atom-surface and the slab-slab configurations. 
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I. INTRODUCTION 

The electromagnetic field interacting with bodies gives 
rise to several effects, among which the radiative heat 
transfer [1] and the Casimir-Lifshitz force [2H1]- These 
two effects are intimately connected, and can be de- 
scribed by a common formalism. If a definite temper- 
ature is assigned to bodies and radiation, two main sit- 
uations are possible. Radiation can be at thermal equi- 
librium with matter everywhere: in this case the heat 
transfer is identically zero, whereas the Casimir-Lifshitz 
force assumes its equilibrium value. If, instead, radia- 
tion is not at thermal equilibrium with matter, both the 
appearance of heat transfer between the bodies and a 
variation of the forces acting on them occur. 

In considering the phenomena of heat transfer and 
Casimir-Lifshitz forces, the bodies shapes, geometric con- 
figurations, and dielectric functions are main issues, to- 
gether with the thermal configuration of the system. A 
general framework used in solving such kind of problems 
is that of macroscopic electrodynamics pW6]. Casimir- 
Lifshitz interaction for systems at thermal equilibrium 
(at T > 0) has been largely studied in the last 60 years 
reaching a consistent and almost complete theoretical 
formulation, as well as experimental observation ^7 . In 
particular, much more recently, the thermal component 
of the force at thermal equilibrium has been measured 
[H] . On the contrary, systems out of thermal equilibrium 
have been much less explored, and mainly simple con- 
figurations and idealized cases (such as infinite bodies) 
have been considered. Nonetheless, systems out of ther- 
mal equilibrium showed remarkable features, already ob- 
ject of theoretical and experimental investigations both 
concerning the force [SHU] and the heat transfer [201 - 
\Tl\ . This is the case of the atom-surface force, where 
the absence of peculiar cancellations among the different 
components of the radiation out of thermal equilibrium 
leads to large quantitative and qualitative modifications 



such as new asymptotic behaviors and possibility of re- 
pulsive interactions [5] . These new features have allowed 
the first experimental observation of thermal effects [12] . 
Recently, motivated by the necessity to develop a more 
complete theory for systems out of thermal equilibrium, 
several studies have been developed [2HHSS]- In partic- 
ular, a general closed-form analytic expression for both 
heat transfer and Casimir-Lifshitz interactions has been 
derived, valid for arbitrary bodies shapes, dielectric func- 
tions, and arbitrary temperatures of the bodies and of the 
environment [28j . To this purpose a scattering-matrix 
approach has been used. This approach has been already 
successfully employed to calculate Casimir-Lifshitz inter- 
actions [16l |28l [Ml |35] : its main advantage with respect 
to the standard Green-function formulation, where the 
electromagnetic problem needs to be solved for the com- 
plete system, consists in requiring the solution for each 
body composing the system independently. 

In this paper we provide a systematic derivation for 
the heat transfer and for the Casimir-Lifshitz force out 
of thermal equilibrium between two bodies at two differ- 
ent temperatures, immersed in and external environmen- 
tal radiation characterized by another temperature |28j . 
Particular attention is devoted to the expression of the 
average values of the electromagnetic field in terms of the 
scattering operators. The general expression is finally an- 
alytically and numerically applied to several simple but 
already interesting configurations, by calculating the heat 
transfer and the force for slab-slab and atom-slab systems 
in an out of equilibrium scenario. The role of finite-size 
effects and of the environmental temperature are shown 
to be qualitatively and quantitatively relevant. 

The organization of the paper is the following. In Sec. 
[TT| we present the physical system, the hypotheses and 
the general theoretical framework. In Sec. |II1| and [rv] we 
evaluate the analytic expression and the fiux of the two 
main ingredients of our calculation: the Maxwell stress 
tensor and the Poynting vector. Section |V] is dedicated 
to the definition of the reflection and transmission op- 
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erators associated to each body. Then, in Sec. |VI| the 
correlators of the fields emitted by the bodies and by the 
environment are expressed as a function of the scatter- 
ing operators. This allows us to calculate, in Sec. |VII[ 
the fluxes of the stress tensor and the Poynting vector in 
any region and then the final expressions of the Casimir- 
Lifshitz force and the heat transfer. This expression is 
first applied to the case of a single body alone out of 
thermal equilibrium in Sec. 



VIII Then, in Sec. IX 



discuss the case of the force acting on a neutral atom in 
front of a planar slab, as well as both the force and the 
heat transfer in a two-slab configuration. We finally give 
some conclusive remarks. 



II. PHYSICAL SYSTEM AND 
ELECTROMAGNETIC FIELD 

Let us start by describing the geometrical configuration 
of our physical system. We are going to deal with two 
bodies, labeled with 1 and 2. From a geometrical point 
of view we will assume that the two bodies are separated 
by a planar surface. This hypothesis is not strictly nec- 
essary in a scattering-matrix approach but it is nonethe- 
less verified in all the relevant experimental configura- 
tions such as for example two parallel planes, a sphere or 
a cylinder in front of a plane, two cylinders, as well as 
an atom in proximity of a planar surface. At the same 
time, this assumption allows us to choose a convenient 
plane-wave decomposition for the electric and magnetic 
fields, leading in this way to quite simple expressions of 
the Casimir-Lifshitz force and the heat transfer. To be 
more specific, the geometry of our system is depicted in 
figure [TJ the bodies 1 and 2 are respectively enclosed in 
the strips zi < z < Z2 and z^ < z < Z4, where Z2 < Z3. 
As a consequence, any plane z = z with Z2 < z < z^ 
separates the two bodies, and three regions A, B and C 
are defined. Our geometrical description coherently in- 




Figure 1: The geometry of the system. The bodies are sepa- 
rated by the strip Z2 < z < Z3. This defines the three regions 
A, B and C. 

eludes as a limiting case the possibility of bodies having 
infinite thickness, taking for example for body 1 (body 
2) the limit zi — >■ —00 (24 — >■ -|-cx)). 



As far as the thermodynamical description of our sys- 
tem is concerned, we will assume that we are able to 
define for body 1 (2) a temperature Ti {T2) and that it is 
in local thermal equilibrium, i.e. the temperature of each 
body is assumed to be constant. Moreover, we assume 
that the two bodies are immersed in a vacuum environ- 
ment (having e = 1) characterized by a third temperature 
T3. We will make the further important assumption that 
the composite system is in a stationary regime, which 
means that the three temperatures remain constant in 
time. 

Let us now first describe the general framework of the 
calculation. We are interested in calculating the force 
F acting on any of the two bodies, as well as the heat 
transfer H on it, defined as the energy it absorbs per 
unit of time. Focusing for example on body 1, these two 
quantities can be expressed under the form of surface 
integrals through a closed surface S enclosing the body 
1 



if = -^(S(R,i))eym 



(1) 



of the quantum symmetrized average of the Maxwell 
stress tensor T (having cartesian components Tly , with 
i,j = x,y,z) and the Poynting vector S. These quanti- 
ties are classically defined in SI units as 

Tij{R, t) = eo [S,(R, t)Ej{R, t) + c'B.iR, t)B^{R, t) 

- ]^(E\R,t) + c^B\R,t)y, 
S(R,t) = eoc^E(R,t) x B(R,t). 



(2) 



and the quantum symmetrized average value (v4i?)sym 
defined as 



IS 



1 



{AB) + (BA) 



(3) 



being (A) an ordinary quantum average value. Before 
working on eq. ([l]) it is evident that a complete descrip- 
tion of the electric and magnetic fields in any region is 
mandatory. 

In our coordinate system the z axis clearly represents a 
privileged direction, being the axis perpendicular to the 
plane separating the two bodies. Inspired by this prop- 
erty, for the mode decomposition of our electromagnetic 
field in any region we replace the common plane-wave 
representation in which a mode of the field is represented 
by the three-dimensional wavevector K — {k^, ky, kz) 
by a description in terms of the transverse wavevector 
k = {kxiky) and the frequency w. The result of this 
choice is that in this approach the z component k^ of the 
wavevector becomes a dependent variable, defined by the 
relation fc^ = — . Since this relation is quadratic in 
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kz, we are obliged to introduce explicitly a variable (f> tak- 
ing the values = ±1 (with shorthand notation = ± in 
the expression of the polarization vectors and field ampli- 
tudes) corresponding to the sign of k^. As a consequence, 
the complete wavevcctor will be noted with 



{kx, ky, (j)kz 



-k2. 



(4) 



For fc < ^ , this relations gives a real value of k^ , and then 
a propagative wave, for which the (f) represents the direc- 
tion of propagation along the z axis. On the contrary, for 
k > kz becomes imaginary and we get an evanescent 
wave whose amplitude depends on z: in this case is 
the direction along with the amplitude of the evanescent 
wave decays. Moreover, we will as usual need an index p 
associated to the polarization, taking the values p = 1, 2 
corresponding to TE and TM modes respectively. Fi- 
nally, in our approach a mode of the field is identified by 
the set of variables (w, k,p, 0). 

The expression of the electric field in any region can 
be first given under the form a frequency decomposition 



E(R,t) = 2Re 



+ 00 



2^ 



exp(— ia;i)E(R, cj) 



(5) 



where a single-frequency component has the following 
mode decomposition 

E(R,c.)=^ j -^cxp(zK*.R)e^(k,w)£;pnk,^) 

(6) 

where from now on the sum on runs over the values 
{-|-, — }, the sum on p over the values {1, 2}. The quan- 
tity (k, w) represents the amplitude of the electric field 
associated to a given mode (w, k,p, 0). For the polariza- 
tion vectors ep(k, a;) appearing in eq. ([6]) we adopt the 
following standard definitions 



The expression of the magnetic field can be easily de- 
duced from Maxwell's equations, and is given by 



where 



exp(iK'^ • R)^p(k, w)i;;J(k,a;) 

(9) 



^;(k,..) = (-i)^4(^)(k,..) 



(10) 



being S{p) the function which switches between the two 
polarization, acting as S{1) — 2 and S{2) — 1. We are 
going to gather the expressions ^ and ([9| of the electric 
and magnetic fields at a given frequency uj in a column 
vector and write 



E(R,tj) 
cB(R,w) 



E 



^exp(zK^.R) 

(27r)2 



(11) 



where in the right-hand side the 1 and the (—1)'' corre- 
spond to the electric and magnetic field respectively and 
we have introduced the vectorial index 



A{p) 



[sIp) 



(12) 



where again the upper (lower) index is associated to the 
electric (magnetic) field. 

Now that the expressions of the electric and magnetic 
fields are explicitly given in terms of a set of field ampli- 
tudes i;^(k, uj) we are ready to work, in the next sections, 
on the explicit expressions of the Maxwell stress tensor 
and Poynting vector. 



e^g(k, cj) = z X k 
etivt(k, w) 



-{-kyic + kxy) 



rj,y^yr.,^) - -e'{^^{k,u) X K"^ = -(-fcz + 0fc^k) 

UJ UJ 



(7) 



where x, y and z are the unit vectors along the directions 
a;, y and z respectively and k = k/fc. The unit vectors 
defined in Q obey the following useful properties 



(-k,^) 



-TE 



(k,^) 



.-0 

^TE 



(k,^) 



(4E(k,' 
(^tm(^' ' 



4E(k,w) 4M(-k,^^) 



4m(]^'^) ^ 

eTM(k,w) kzi 



(k,c.) 



(8) 



III. THE MAXWELL STRESS TENSOR 



A. General expression of the tensor 

From the definition of the stress tensor ^ it is clear 
that we need to calculate explicitly the quantities EiEj 
and BiBj for i, j = x,y, z and as a further step their sym- 
metrized quantum average. Using the compact vectorial 
notation introduced in eq. (11), we have 
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where the fields imphcitly depend on (R, t) and Z gener- 
icahy gathers all the terms proportional to EE or E'''E^ 
whose average quantum value is zero. Having in mind the 
transition to the quantum symmetrized average values, 
we now give the following definition for the commutators 
of the field amplitudes 



{E^ik,u;)Ef,Hk',u;%,^, = ^{E^ik,u;)Ef,\k' ,u;') 

+ Ef,H^',u;')E^{Ku^)) 

= 2TrS{Lj-Lj'){p,k\C^^'\p\k') 

(14) 

where we stress the fact that in general two modes of the 
field propagating in opposite directions do not necessarily 
commute. Moreover we have explicitly inserted the con- 
servation of frequency: since in our system no dynamics 
is considered, the field amplitudes at different frequencies 
necessarily commute as a consequence of time invariance. 
In virtue of this conservation any correlation function of 



the electromagnetic field analogous to ( 14 1 is function of 
a single frequency w (and not of both w and w'), which 
will not be explicitly written from now on. Moreover, 
in eq. ( 14 ) we have expressed the correlators as matrix 



elements of a matrix C'*"'' : in our notation, the matrices 
are defined on the space {p, k) being p = 1,2 and k G 
and thus the product of two matrices A and B is given 
by 



{p,k\AB\p' X) 



E 

p" 



d'k 



(2^) 

X (p",k"|6|p',k' 



^{p,k\A\p",k" 



(15) 



where the matrices A and B are both calculated at the 
same fixed frequency lo. 



Using eq. ( 13 ) and the definition ( 14 1 we have 



EiEj 



EE 

tfnj)' pp' 



£k 



(fk! 



+°° diO 
2tj: 



(— 1)^+^"' 



exp[i(k — k') • r] exp[i((/)fcz — (/)'/c*)z] 



(p, k|C^^' b', k') { (k, c)) ^ (4;;,) (k', u:)) 



(16) 
defined in 

( 14 1 and appearing in ( 16 ) depends on the region (A, B 



It is important to note that the correlator 



or C) in which the average is calculated: this information 
is contained in the z dependence of the fields in the left- 
hand side of eq. (16), as well as in the z coordinate 



explicitly present in its right-hand side. As we will see in 
section [ill B[ we only need to evaluate the symmetrized 
average value of the flux of the iz components {i = x,y, z) 
of the electromagnetic stress tensor. These quantities can 
be calculated using eqs. ^ and (16 1. In the case of Tzz 
we have 



(2^ I (2^ -exp[z(k-k')-r]exp[z(#,-0'fc;)z](p,k|C^^b',k') 



x{(^^f{k,u))^[ef,*{k\L.))^-[i^^{k,c.))^{e^^^^^ 

+ (- ir ''I (4(p) (k, -)) ^ (k^ -)) ^ - (4(p) (k, -)) ^ (4i 



(17) 
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whereas for T^z (being m = x,y) a straightforward calculation gives 



(2^ / (2^ / 2^ ~ ^'^ ■ ^^pW<^^^ - k|C*^ b', k') 



{ ™ (^t' * (k'. -)) ^ + (-1)^"^' (4(p) (k, -)) ^ (4;;, (k', -) 



+ rn ^ z 



(18) 



where the symbol m <=^ z represents two more terms ob- 
tained by interchanging in the first two m with z. Let 
us make some comments on these expressions. First we 
note that, even if the stress tensor is in general a func- 
tion of time, this is not true for its quantum average 
value. Besides, this is a function of the correlators of the 
electric-field amplitudes (propagating either in the same 
or in opposite directions) and of the components of the 
polarization unit vectors. As we will see in the next sec- 
tion, the components (17) and (18) are the only ones we 



need to calculate all the components of the force acting 
on the two bodies. 



B. Flux of the stress tensor in terms of field 
correlators 

As anticipated, the electromagnetic force acting on a 
given body inside a volume region can be calculated by 
taking the flux of the stress tensor through a closed sur- 
face enclosing this volume. As a consequence, the m 
component of the force on a given body is given by the 
flux 



F„ 



{^mj ) syin 



(19) 



where the summation over repeated indices is assumed 
and E represents any closed 602; entirely enclosing the 
body. Let us choose for example the box depicted in fig. 
[2] i.e. a parallelepiped having one side of length D and 
as a base orthogonal to the z axis a square of side L. 
According to the definition of T^j, the m component of 
the force (m = x, y, z) is given in this case by the flux of 
Tmz through the two surfaces orthogonal to the z axis, 
plus the fluxes of T^x and T,„j, through the surfaces of 
the parallelepiped orthogonal to the x and y axes respec- 
tively. Taking now the limit L -l-oo we see that the 
surface of the two bases orthogonal to the z axis diverges 
more rapidly (like I?) than the other four surfaces (like 
L). As a consequence, we deduce that in order to cal- 
culate Fra one simply needs to calculate the flux of T^z 
on the surface (which has now become a plane) in region 
A and subtract this result from the flux of Tmz through 
the plane in region B. Moreover, due to the arbitrariness 
of the box, these two fluxes must not depend on the z 
coordinates of the respective planes, even if in general 
the stress tensor depends on z. 

From this discussion we conclude that we need the flux 
of Tmz through a plane z = z. We will then have z < z\ 



D 



B 



Z2 



Figure 2: The box chosen for the calculation of the force acting 
on body 1. 



for region A, Z2 < 2^ < -2:3 for region B and z > 2:4 for 
region C. Integrating eq. (17) on the plane z = z and 
noticing that this gives a Dirac delta (27r)^(5(k — k'), we 
get the flux of T^^ expressed as a function of the field 
correlators 



^z(z)= f (fv{Tz, 

'J Z — Z 



eoEX! / / !!r exp[i(#3 - (/)'fc*)z] 

■j>i pp' 



(2^)2 7„ 27r 



^{i't)^mz-('t)Mi.-i't)M'). 

-(^i(p)),«;'))j}(^'kic^^'b',k) 



(20) 



where all the polarization unit vectors are calculated in 
(k, Ld). For the other two components of the stress tensor 
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(m = x, y) we have, integrating eq. ( 18 1, 

J Z — Z 



EE 

pp' 



I sym 

2^ 



(2^) 



e:xjp[i{<f>kz — (j)'k*^)z\ 







+ m ^ z 



}(p,k|C**'b', 



k). 



(21) 



From the definitions ([t]) of the polarization unit vectors 
we deduce that 

+(-i)^"1fe)),(4;;o).-(4(.)),(4;;.). 



(22) 



and 



We now observe that 



(23) 




2(/)Re 
2i(blm k. 



(24) 



For both expressions, in the former case (0 = (j)') only 
the contribution coming from propagative waves plays a 
role, whilst in the latter only evanescent waves are rele- 
vant. At the same time, for both expressions, and both 
for <j) — (j)' and (j) ^ (j)' , the exponential term contain- 
ing z disappears, as expected: whereas the stress tensor 



depends on z, this is not the case for its flux calculated 
on a plane having an arbitrary position z = provided 
that z — z remains in a given region (A, B or C). Finally, 
the flux of Tmz on the plane z — z can be cast for any 
m = x,y, z in the form 



duj 
2tt 



E 



ck 



, Jo 



2eoc2fc, 



(p,k|C^^ |p,k) X 



m — z 



duj 
2tt 



(25) 



This equation represents the main result of this section, 
describing the flux of the component Tmz of the stress 
tensor as a function of the field correlators. As we have 
shown before, this quantity, together with the explicit 
knowledge of the correlators in any region of our system 
which will be discussed in Sec. |VI[ is sufficient to deduce 
any component of the force acting on the two bodies. 
From eq. ( 25 ) we deduce that the flux is written as the 



sum of two separate contributions, coming from the prop- 
agative and evanescent sectors respectively: the former 
depends on the correlators between the field propagat- 
ing in a direction (p and itself, whilst the latter implies 
the correlators of counterpropagating fields. We also re- 
mark that the quantity C"^* is the only term in eq. ( 25 ) 



depending on the region in which the flux is calculated, 
which means on the position of z. Finally we observe 
that the result deduced for the flux of Tzz {'m = z in eq. 



(25)) coincides with the expression obtained in jT5| . 



IV. THE POYNTING VECTOR 

A. General expression of the vector 

Let us now focus our attention on the Poynting vector 
defined in eq. In order to evaluate its quantum 

symmetrized average we first need to work out the generic 
field product EiBj. Using the same conventions of the 
last section we obtain 



EE 

4><. 

|exp[i( 



' pp' 



d^k f d^k' 



duj 



(27r)2 J (27r)2 2it Jq 2tt 



K'^'*).R]expH(^-a.')i](-l)^' f ^^f' ""^^^''l'"'' ""'^ 1 fe(k,a.)) (4*„{k',^ 

) i yi y J iy ) yEy{k',u')E^{k,iu)i y''^ ' VA ^^(p')^ ' 

E^^k,oj)Ef,{k',u;')' 
E'^:ik\u')Ep{k,u:), 



^Ey{k',u;')E^{k,u;) 

+ exp[.(-K^* + K>') . R] exp[.(. - .'M-lf fS^' ^.^^'^^klll (s^* (k,..))^^^^^''-')) ,} + 



Z 

(26) 
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from which we immediately get 

d^k f (fk' r+°° duj 



{cE,B 



EE 



(t>(t>' PP' ^ / ^ V J JO 

X exp[i(k — k') • r] exp[i((/)/c^ — 

X {p, k\Cf"f' \p', k') { (-1)^'' [4{k, c)) ^ (4;;,) (k', 

+ (-ir(4(^,(k,.)).(.t:*(k',.))j 



(27) 



with the same conventions of eq. ( 13 ) 



B. Flux of the Poynting vector in terms of field 
correlators 



We now observe that, in virtue of the same discussion 
about the closed surface S given in the last section, we 
only need the flux of the z component of the Poynting 
vector in order to evaluate the heat flux on one of the two 
bodies. For the flux of on the plane z = z we have 



d^k 



= ^ocEE 



(2^)^ Jo 



+00 



duj 
2^ 



[x^y]}{p,k\C'^'^'\p',k) 



(28) 



where we stress the fact that the two terms obtained by 
interchanging x and y must in this case be changed in 
sign. We then obtain 



(-i)^'Gt)i4'(;o),+(-i)''(4(.) 

= Spp' — {(l)kz + (l)'kl) 



(29) 



and finally cast the expression of the flux of the z com- 
ponent of the Poynting vector under the form 




2eoc2#2 



(p,k|C*^'|p,k) 



(30) 



This expression is the analogue for the Poynting vector of 
eq. (25). It constitutes, together with the expression of 
the matrix elements of C^^ at any frequency, the main 
ingredient in the calculation of the heat transfer on the 
two bodies. 



SCATTERING FORMALISM: REFLECTION 
AND TRANSMISSION OPERATORS 



The quantities C"^* appearing in eqs. (25) and (30) 



and defined in eq. (14) are the correlators of the total 



fields in each region. In our problem these fields result 
from the ones emitted by the two bodies and the envi- 
ronmental field, as well as from all the possible scatter- 
ing processes undergone in presence of the two bodies 1 
and 2. We then need to introduce a set of operators de- 
scribing the scattering produced in presence of a single 
arbitrary body. Let us suppose to have a body located 
in the region zi < z < Z2. Let us further assume that 
an external field is impinging on our body, either from 
its left or from its right side. This field will be scattered 
upon the body, producing in this way new components 
of the field on both sides of the body. In particular, the 
field coming from the left (right) will produce a reflected 
field propagating toward the left (right) on the left (right) 
side, and a transmitted one propagating toward the right 
(left) on the right (left) side. The two possible configu- 
rations, depending on the direction of propagation of the 
incoming field, are represented in Figure[3j The reflection 




T- 



n 



+ 



7(tr) 



Figure 3: The definition of the reflection and transmission 
matrices. 



and transmission matrices TZ^ and are the operators 
linking each mode of the outgoing fields to the incoming 
ones. 

In particular, considering the case of a field coming 
from the left side, the incoming field 



E('")+(R,i) = 2Re 



V- /■+°° doj f d'^k , . 



exp[-^a;^]e+(k,CJ)^;('"^+(k,w) 



(31) 



will result in a reflected (on the left) and a transmitted 
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(on the right) field defined by 



VI. FIELD CORRELATORS 



(R,0 = 2Re 



exp[-iwt]ep(k,t^)4'-^)-(k,w) 



(R,t) = 2Re 



E 



27r 



(2^)2 



(32) 



exp[iK+ • R] 



exp[-zwi]e+ (k, w)S^*^)+ (k, cj) 



(33) 



As remarked before, since the scattering process is sta- 
tionary, the frequency is conserved. We are thus able to 
define the operators TZ~ and through the following 
relations involving the amplitudes defined in eqs. (31 1, 
(1321) and (ISSl) 



V 



-(k,c.) 



E 

p' 



fk!_ 
(2^ 



(p,k|7e-b',k')4r+(k',^) 

(34) 



i?f)+(k,c.) = 



p' 



(p,k|r+b',k')4!"'+(k',u;) 

(35) 

connecting each mode of the outgoing fields to all the 
modes of the incoming one at the same frequency lo. In 
analogy with the field correlators, we are going to drop, 
for the sake of simplicity, the dependence on the fre- 
quency LO in the reflection and transmission operators. 
A perfectly analogous procedure leads to the definition 
of the scattering operators TlT and ■ 

Using the formalism we have just introduced, the case 
of the absence of a given body is obtained, as far as its 
scattering operators are concerned, by imposing 



= 



= 1. 



As we will see later, it is convenient to introduce a mod- 
ified transmission operator which, in analogy with the 
reflection operator, goes to zero as well in absence of the 
body. We thus define 



7-0 = 1 + 7-0 



(37) 



writing the transmission operator as the sum of the iden- 
tity, describing the incoming field propagating unmodi- 
fied on the other side of the body, and of a new operator 
T'^ accounting only for the scattered part of the field. In 
the limit of the absence of the body, we have as desired 
T"^ = 0. 



In order to proceed further and to calculate the force 
and heat transfer on body 1, we need an expression for 
the correlators (p, klC"^"^ |p',k'), defined in eq. (14), in 
each region of our system. These correlators will be ex- 
pressed as a function of the correlators of the field emit- 
ted by each body by means of the scattering operators 
introduced in Sec. [V] In this section we will first address 
the characterization of the single-body field correlators 
(with this expression we refer both to the bodies 1 and 2 
and to the environment), and immediately after express 
the total-field correlators as a function of these quanti- 
ties through the scattering operators. We will note with 
i?p^-''^(k, w) each mode of the total field propagating in 
direction (j) in the region ^ — A,B, C , as shown in figure 
[1] In order to calculate the fluxes ( 25 1 and ( 30 1 in the 



three regions, we need to know the expression for the the 
correlators C;^*^ defined by 

(4^'^k,w)47)^'^(k',L.'))sy,n = 2t:5{u:~J) 

X (p,k|Cf b',k'). 



(38) 

For a system as the one in figure 1 at thermal equilibrium 
at temperature T with the environment, the correlators 
of the total electromagnetic field outside the body follow 
from the fluctuation-dissipation theorem [3 6) 

(i?f°')(R,w)£;f°*^^(R',c.'))sym = 2i,5{lo-lo')-N{lo,T) 



X ImG'y (R,R',w) 



(39) 



being 



N{u:,T) = ^coih[-^)=hu^[^+n{u:^T)\ (40) 



with 



n(w,T) 



e^BT _ 1 



(41) 



(36) In eq. (39) i and j refer to the cartesian components 



of the field and Gy(R, R',a;) is the ij component of the 
Green function of the system, solution of the differential 
equation (see also appendix O 



Vr X Vr - ^e(tj,R) 



(R, R', w) = — ^ M(R - R') 

(42) 

being I the identity dyad and e(w, R) the dielectric func- 
tion of the medium. The property ( 39 1 does not hold in 



the case of a general nonequilibrium configuration. In our 
particular system we have assumed that for each body a 
local temperature can be defined, and remains constant 
in time. This assumption reasonably leads to the hy- 
pothesis that the part of the total field emitted by each 
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body is the same it would be if the body was at thermal 
equilibrium with the environment at its own tempera- 
ture. In other words the emission process is not con- 
siderably influenced by the modification of the external 
radiation impinging on the body. This hypothesis implies 
that the correlators of the field emitted by each body can 
still be deduced using the fluctuation-dissipation theo- 
rem eq. (39) at its local temperature. We note that 



the limits of validity of this hypothesis, already used in 
[U O [T3l [Ml |37] , require further experimental and theo- 
retical investigations. 

We are now ready to write down the expressions of the 
correlators of the environmental field as well as the fields 
emitted by the bodies at local thermal equilibrium. 



A. Correlators of radiating bodies and environment 

1. Environmental field 

The correlators of the environmental radiation in equi- 
librium at temperature T3 are well known: they are given, 
for e {+,-}, by 



S^,4>' w^N{uj, T3) Re(^) Spp, i2n f6iuj - uj')Sik - k') 



2enc 



5^.^'27r<5(w-w')(P,k|C(3)|p',k'). 



(43) 



In this expression we have defined the matrix C^'^-' and 
introduced the notation, valid for any integer n, 

(p,k|7'iP"/™)b',k') =A:;'(p,k|n(P^/™)|p',k') (44) 

n(p*) (II^''*)) being the projector on the propagative 
(evanescent) sector. We remark that the operators 
.p^pw/cw) -|-|-(p.j,,yp^-) |j^gpgj_^|j^ implicitly on the frequency 

LO. 



2. Field emitted by eacti body 

We now calculate the correlators of the field emitted 
by body i (i = 1,2) at temperature Ti by ignoring the 
presence of the other body and assuming thermal equilib- 
rium at temperature Tj . The main point of this deriva- 
tion is the connection between the Green function and 
the scattering operators: the details of this calculation 
are presented in appendices [C] and |D] The result for the 
field correlators is, for two modes of the field propagating 



in the same direction, 

(£;«^k,..)£;«^^(k',..')).ym = . 

ZtQC 

X 2ndiLj - Lo'){p, k| (t'It' - 7^(')'^7'i7^7^(^)^^ 

(45) 

being T?,'*-'"^ and T^'-''^ respectively the reflection and 
transmission operators associated to the side (j) of body 
i defined in section |Vj For fields propagating in opposite 
directions (0 ^ cj)') we have 



(i?«nk,c.)i?«*'^(k',..')) 



X 2TrS{uj — uj' ){p, k| (^-7^(*)'^7?^P|")r^*^"*'^ 

_ -J-{i)4"p(P^)'Jl(i)<i>'->i _|_ 7-(»)<!ip(cw) _ ^(ow)^(i)0'|^ |p',k'). 

(46) 

In analogy with the previous definitions, the correlators 
of the field produced by the body i will be gathered in 
the matrix (7^*^'^'^ , defined by the relation 

(£;«nk,^)4:^*'^(k',c.')).y„. = 27r<5(u;-a.') ^^^^ 

X (p,k|C«^^'b',k'). 

We are now ready to characterize the total field in each 
region by means of the scattering operators, and then to 
deduce its correlators using the results just obtained in 
this section. 

B. Correlators in region B 

In order to build up the field in the region B between 
the two bodies the ingredients we need are the amplitudes 



i;^^^+(k,w), E^-''~(k,u) and £;^''^=^(k,a;). In the region 
B of figure [l] the field propagates in both directions: its 
amplitudes will be simply noted with i?p^'''^(k, w). Gath- 
ering all the modes Ep^^'^{'k,uj) in the symbol i?*^^^"^, the 
amplitudes can be expressed as the solutions of the sys- 
tem of equations 



,(2) 



(3)±/ 



^(B)+ ^ 



7-(i)+£;(3)H 

-7-(2)-£'(3)- 



.7^(l)+s(B) 

.7^(2)-^(B) 



(48) 



where all the operators and field amplitudes are calcu- 
lated at a given frequency w, not explicitly indicated, 
and the products between scattering operators and fields 
are to be considered as matrix- vector products. As an 
intermediate step, we have 



'^(B) + 



= ^7(2l)7^(2)-^;(l)+ + [/(2i)7^(2)-7-(i)H 

+ C/(21)7-(2)-£;(3)- u(21) £1(2)- 



i5(3) + 



(49) 
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where we have introduced the operators 



of the scattering operators of the two bodies 



c/(i2) = (l_7^(l)+7^(2)-)-l 
c/(2i) = (l_7^(2)-7^(l)+)-l 



(50) 



describing the series of intracavity (between the two bod- 
ies) reflections produced by the single-body operators 
7?.*^"'^^"'" and 72.'^^-'". From the definition 



t/(i2) ^ (1 _ 7^(i)+7e(2) 



+00 



5](7^(l)+7^ 



(2)-Nn 



(51) 



)^(k,w)4^')^'^(k',L.'))syn. = 2nSiu; - u;') 



(B)0't 



(56) 



and its analogous counterpart for C/^^^^ we easily deduce 
the following useful properties 



7^(i)+[/(2i) ^ [/(l2)^(l)^ 
and 



7^(2)-[/(12) ^ jj(21)-^(2)- 
(52) 



^(l) + f^(21)^(2)- ^ ^(12) _ ^ 
^(2)-^(12)^(l)+ ^ ^(21) _ 



(53) 



These relations allow us to obtain the following final ex- 
pression of the field propagating in both directions in 
region B as a function of the fields emitted by the bodies 
and the environment 

'^(B)+ ^ jj(12) ^(l)-l- _,_ 7^(1)+(7(21) ^(2)- 

+7^(l) + (7(21)7-(2)-£;(3)- ^ (7(12)-y-(l)-(-£;(3) + 
^(B)- ^ ^(2)-^(12)^(l)+ ^ jj{21) ^(2)- 

+ [/(21)7-(2)-^(3)- _,_ 'j^(2)-jj(12)q-{l) + £1(3) + 

(54) 

We remark here that by taking E'^^^^ = or T"*-^-*^ — 
7^(2)- = in eq. (54) we go back to cqs. (18) and 
(19) of [H], where for both bodies infinite thickness was 
assumed. Since the fields E^^^^^ and E^^'>^ clearly do not 
participate in the expression of the total field between 
the two bodies, the expression of £'(^)'^ (for (/) — +,—) 
can be cast without loss of generality in the form 

£(B)± ^^(B)±£,(l)+^^(B)±^(2)-^ J2 5(B)±£(3)a 

(55) 



which in this case gives, by comparison with eq. ( 54 ) , 



^(B)+ ^ c;(12) 



^(B)+ ^ ^(l) + f;(21) 



^(8)+ ^ .^(l) + j^(21)-y-(2)- ^(B)+ ^ JJ{12)^{1) + 
^(B)- ^^(2)^^(12) 



B 



(B)- 



B 



(B)- 
2 

(B)- 



= 7^(2)-[/(l2)7-(l)+, 



Using eq. ( 55 1 the correlators in region B can be ex- 
pressed as a function of the correlators C'*-'"^'^ (for i = 
1,2) and C(3) given by eqs. Q, (|45]) and (|46]). We fi- 
nally obtain the expression of the matrix Cg"^ in terms 



C. Correlators in regions A 

The complete knowledge of the properties of the field 
in region B is not sufficient, in general, to deduce the 
force and the heat transfer associated to any of the to 
bodies. Focusing our attention on body 1, for example, 
we also need to characterize the field in the region on its 
left side, namely region A. The field propagating 
toward the right in this region is obviously only the en- 
vironment field propagating in the same direction 
On the contrary, as far as the field is concerned, 

it will also include, assuming a finite thickness for body 

1, components from the fields produced by bodies 1 and 

2, as well as from the environment field i?^'^)". The total 
field in region A is then entirely described by the system 
of equations 



£(A)- ^ £(!)- ^ 7^(l)-s(3)+ + 7-(i)-£;(B)- 



(57) 



which using the result (54) for the intracavity field be- 
comes 

-£;(A)+ ^ ^{3) + 
£;(A)- ^ £;(!)- +7-(l)-7^(2)-[/(12)£;(l) + 

+7-(l)-[/(21)£;(2)- j(l)-u(2Dq-(2)- -^(3)- 

+(Tzw- + r(l)-7^(2)-{/(l2)r(l)+)£;(3)+ 



Using the general decomposition 



(58) 



^ ^ 4^)ii;(i)"+4^)£;(2)-+ ^ B(A)£(3)a 

Q— +,— a— 4-,— 

(59) 

we have 



^(A) ^^(l)-^(2)-^(12) 



A[t^ = 1 



(A) 



A) 



7^(1)- +7-(i)-7e(2)-[/(i2)7-(i)H 
7-(i)-(7(2i)7-(2)-^ 



We are now ready to give the final expression of the cor- 
relators C^"^ of the total field in region A in terms of the 
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scattering matrices, which reads 



bodies and the environment ( 43 ) , ( 45 ) , ( 46 ) 



C++ = C(3) 

c+- - 

(60) 



Due to the geometry of our system, the correlators of 
the field in region C can be obtained from the ones given 
here for region A performing the interchanges A ^ C, 
1 2 and + ^ — . This holds for all the other quantities 
we are going to calculate in region A. 



VII. FINAL EXPRESSIONS FOR THE 
CASIMIR-LIFSHITZ FORCE AND HEAT 
TRANSFER OUT OF THERMAL EQUILIBRIUM 



Casimir-Lifshitz force 



We are now going to calculate the flux of the stress 
tensor, in order to deduce the expression of the force. For 
simplicity, we will focus on the z component of the force 
acting on body 1, using as a consequence the expression 
(251 with TO = in regions B and A. The calculation 



of the other components of the force follows the same 
scheme we are going to present in the following. Let us 
now turn, then, to the evaluation of the fluxes of T^z in 
regions B (Sec. |VII A l\ and A (Sec. |VIIA2[): wc will 



assume that for both surfaces the vector orthogonal to 
the surface is oriented toward the right, i.e. the positive 
direction of z axis. These two results will provide the 
final expression of the Casimir-Lifshitz force acting on 
body 1, deduced in Sec. |VIIA3[ 



1. Flux in region B 



Using the relation ([25| (with to = z and Z2 < z < z^) 
and ( 56 ) , we are able to express the flux in region B as 



X (p,k|{n(p-)(^p)+c(i)++^(B)+t 

Q = + ,- 

+ n(™) (s(B) + c(3)i3(B)-t + 5(B)- c(3)5(B)+t^ ] I k) . 

(61) 

Defining the trace operator for a frequency-dependent 
operator A as 



we can write the flux in region B under the form 



(62) 



P. 



(pw) 



$(^H7^i,r2,r3) = -2eoc2Tr|^ 

+ Af )+C(2)-Af + Af '-C(2)-Af )-^) 
+ Af )+C7(2)-Af + Af )-C(2)-Af )+^) 



a=+,- 



(63) 



In both expressions it is clear that we have three sep- 
arate contributions associated to body 1, body 2 and 
environment respectively. Using the fact that the trace 
is invariant under cyclic permutations, we have 



$(^)(ri,T2,T3 



Tr|l 



a function of the correlators of the field emitted by the 



X N{uj, Ti)J{n^^^+, 7^(2)-) + N{u, Ts) J(7^(2)-, 7^(l) 

+ (7V(w,r3) - A^(t^,Tl))i^(7^(l)+,7^(2)-,r(l'+) 
+ (7V(c^, Tg) - Nicj, T2)) if (7e(2)- , 7^(l)+, r(2)-)] } 



(64) 
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where 

J(7^(l)+,7^(2)-) = ;7(i2) (^7,(pw) _ ^(i)+^(pw)^(i)+t 

J/(7^(l)+^7^(2)-^7-(l) + ) ^ ^(12)-y-(l)+ppw^(l)+tf^(12)t 

(65) 

Note that when calculating J{71^^^~ , T?.*^^^"*") one also need 
to change C/^^^^ into C/^^^l 

Eq. ( 64 1 can be cast in the form 



$(^)(ri,r2,T3) = 

+ A$(B)(ri,T2,T3) 

where 

= -2ReTr|^7V(c^,T) [[/(l2)7^(l)+7^(2)- 



(66) 



(67) 



and 

A$f)(Ti,r2,r3) = -s 



X Tr 



"42 

2 



(J(7^(l)+,7^(2)-) - J(7^(2)-,7^«+)) 



+ n3li^(7^(l)+,7^(2)-,r(l)+) 
+ n32i^(7^(2)-,7^(l)+,r(2)-)" 

where we have defined, for i,j = 1, 2,3, 
Uij = n{Lo,Ti) - n{uj,Tj). 



(68) 



(69) 



At this point, a remark is important about the expres- 
sion (67 1, giving the equilibrium part of the flux in re- 



gion B of T^^. We have to observe that, since the operator 
is diagonal in the (k,p) basis, its trace defined as in 



,(pw) 



(62 1 is divergent. Moreover, this term is independent of 



the bodies under scrutiny. Nevertheless, we have to keep 
in mind that the flux in region B does not have a direct 
physical meaning, since we still have to subtract from it 
the flux in region A in order to obtain the value of the 
force. As we will see in Sees. VII A3 and VII C[ the di- 
vergences present in the individual fluxes are completely 

regularized when taking the difference ^i^'' — 



2. Flux in region A 
In order to calculate the flux of T^z in A we have to 



use eq. (251 with m — z and z < zi. Moreover, the 



correlators in region A are given in eq. ( 60 ) . After al 



gebraic manipulations analogous to the ones used in the 
last section we have, for the flux in region A, 

<i>(A)(ri,T2,T3) = ~2eoc2Tr[^7'^P*)(c(3) 

+ s(A)c.(3)5(A)ty ^ 

(70) 

This can be cast in the form 

(71) 



where 

A$('^)(ri,T2,T3) = -HTt 
+ ?^32'P2^'"''T^^^"C/'^^^T^^^"7'iP"V^^'"'fC/'^^^^T*^^^''' 

+ "21 (t'IT' + 7^(2)-7'i7' - 7'i7)7^(2)-t 

_ 7^(2)-7?^P«')7e,(2)-t^ [/(21)t7-(l)-tp^Pw)7-(l)-[/(21) 

+ -p^p™)7-(i)-7^(2)-[/(i2)7-(i)+-p(F")7^(i)-t 

_|_-p^P*)7-(l)-7^(2)-[/(12)-y-(l) + -p(pw)-y-(l)-|-t 

X [/(I2)t7^(2)-t7-(i)-ty . 

(72) 

We repeat here that the flux in region C (necessary for the 
calculation of the force acting on body 2) can be obtained 
from eqs. (71 1 and (72) by performing the interchanges 
A ^ C, 1 ^ 2 and -h ^ -. 



3. Casimir-Lifshitz force acting on body 1 

We now have all the ingredients to give the z compo- 
nent of the force acting on body 1. From the definition 
of the stress tensor we have 

Fi, = <i>''P\Ti,T2,T3) - $(A)(ri,T2,T3) (73) 



where the two fluxes are given by eqs. (|66|) and ([71| 
Gathering all the results obtained in the previous sec- 
tions, the complete expression of the force reads 

TTi(oq)/'T^ \ , Ci(eq)/rTi \ 
p /rp rr. X (J-l)+-fz [J- 2 ) 

-riz(Jl, J2, -tSj — 7, 



(74) 



■AFUTuT2,T3). 
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In this expression the resuh is written as a sum of two 
terms. The first contribution is the average, at the tem- 
peratures Ti and T2 of the two bodies, of the equihbrium 
force 



where 



Fi'='5)(T) = -2ReTr|— 7V(a;,r) 

} 



C/(i2)7^(i)+7^(2) 



^(21)^(2)-^(l) 



(75) 



which contains both the zero-temperature term and the 
thermal correction. This result for the equilibrium force 
was already obtained by different authors in the frame- 
work of scattering-matrix theory [MJ |3S] . As remarked 
in [IB] , the eq. ( 75 ) gives a finite result for any choice of 
temperature and material properties for the two bodies. 
Moreover, the equilibrium force (75 1 shows the important 



property of depending only on the intracavity reflection 
operators 72.'^-*+ and TZ^^^^^ , i.e. the operators describing 
the reflection produced by each body on the side of the 
other one. 



The second term in ( 74 ) is the non-equilibrium contri- 
bution, given by 

AFi3(Ti, T2, Tg) = A$P) (Ti, T2, Tg) - A^^^) (Ti, T2, T3) 
+ ;iTr(n327'lP"') 

(76) 

where the two fluxes A'&i^'' and A^i^^ are explicitly 



given by eqs. (68) and (72 1 respectively. The non- 



equilibrium contribution manifestly satisfies the condi- 
tion 



AFi,(T,T,T) =0. 



(77) 



Differently from the equilibrium force ( 75 ) , the nonequi 



librium contribution (761 still contains terms which are 



individually formally divergent. In Sec. VII C where a 
unified expression for the Casimir-Lifshitz force and the 
heat transfer will be provided, we will see that this can be 
manipulated so that all these divergent terms disappear. 



B. Heat transfer 

In order to obtain the expression of the heat transfer 
on body 1 we have to follow the same steps we used in the 
case of the force. We first need the fluxes of the Poynting 
vector in regions B and A. Their difference will provide 
us the energy absorbed per unit of time by body 1. The 
flux in region B can be obtained by combining eq. ( |30[ ) 
with the correlators in region B given by eq. ( 56 ) . The 
result can be cast under the form 



^) (Ti, T2, T3) = hTT{uj [n3l/:(7^(l)+, 7^(2)-, r(i)+) 
n23£(7^(2)-,7^w+,r(2)-) + n,2j{n^'^+,n^^'>-)] } 



)+t 



:7(7^(l)+,7^(2)-) = [/(12) (^^^m _ 7^(i)+p(pw)^(i 

X (p(P'') _ 7^(2)-t75(p-)7^(2)- 

^^(2)_tp(™) _^(cw)^(2)-^ 
£(7^(1)+^ 7^(2)-^ 7-(i) + ) = ^(12)7-(l) + -pPYX(l)+t[/(12)t 

X (ri"""^ - 7^(2)-t7?(p*)7^(2)- 

(79) 

An analogous calculation leads us to the following ex- 
pression of the flux in region A 

+ n2i (t'IT^ + 7^(2)-75i7) - 75(<=-)7e(2)-t 

_ 7^(2)-p(P|")7^(2)-t^ [;(21)t-y-(l)-t.p(pw)^(l)- j^(21) 

+ (7'(p")7^(l)-7'i7'r(^)+^f/(^')^7^(2)-tr(l)-^ 

_j_ -p(pw)^(l)-.^(2)- [;(12)-y-(l) + -p(pw).^(l)_t 
_^7r,(pw)-y-(l)--^(2)-^(12)-y-(l) + p(pw) 

X r^i)+t;7(i2)t7^(2)-t7-(i)-ty ^ 

(80) 

The total heat flux on body 1 is finally given by the 
difference of the two contributions 

i7(ri,T2,T3) = (^(AHri,r2,T3)-(^(B)(Ti,r2,r3). (8i) 



Since the fiuxes ( 78 ) and ( 80 1 in regions B and A rcspco 



tively are zero for Ti = T2 — T^, the heat flux (81 1 on 
body 1 satisfies the evident property 



H{T,T,T) = 
for any temperature T > 0. 



(82) 



(78) 



C. Unified expression for force and heat transfer 

We are now ready to give the main result of the paper, 
namely the analytic explicit expressions of the Casimir- 
Lifshitz force and heat transfer on the body 1. These 
expressions are valid for any choice of the shape and di- 
electric properties of the two bodies. We are going to 
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give the foUowing definitions 

Fi. (Ti , , ) = ' 2 + A2 (Ti , T2 , T3 ) 

(83) 

H{Ti,T2,T3) = Ai(ri,T2,r3) (84) 

we can coUcct and to give a unified expression for the 
nonequihbrium contribution A2 to the force and the heat 
transfer Ai, both relative to the body 1. Before provid- 
ing the explicit analytic expression of A^ for m = 1, 2 
we recall that the fluxes (71) and (66 1 of T^z in regions 



A and B respectively contain individual divergent terms. 
The same property holds for the heat transfer, as it is 
evident from example from the flux ( 80 ) of Sz in region 
A, containing as a first term the trace of n^P*^ We are 
going to show that the nonequihbrium force and the heat 
transfer are indeed convergent for any choice of the two 
bodies. To this aim a fundamental intermediate step is 
the identification of the individual divergent terms in the 
expressions of the fluxes of Tzz and Sz- We first observe 
that all the terms which do not contain any reflection or 



transmission operator, as the ones we have already dis- 
cussed, are indeed divergent. This in not the case, on the 
contrary, for the terms proportional to at least one reflec- 
tion operator, since these ones tend to zero in absence of 
the objects. As far as the transmission operators are con- 
cerned, we have then to express each T operator as 1+T: 
in analogy with the reflection operators, T tends to zero 
in absence of the bodies. Finally, considering the terms 
containing only the operators C/'-^^-', U^"^^^ and projection 
operators, it is sufficient to use the relations (53) in order 



to write each of them as a sum of a divergent one, which 
is independent on the scattering operators, and another 
one proportional to the reflection matrices. 



By following the procedure we have just described it 
can be shown that all the divergent terms exactly cancel 
each other. All the remaining terms are proportional to 
either a reflection or a modified transmission operator, 
as explicitly shown in eqs. (85)-(88) below. We are now 



ready to give the final analytic unified expression for the 
Casimir-Lifshitz force and heat transfer on body 1. This 
reads 



A™(ri,r2,T3) = {-ir+'hTr [(c/(2i)t(2.g,„(r(i)-) - U{n^'^+))u^''^ + u„,) (pi^) + /-l(7^(2)- 

+ (-l)™(c/(l2)/_l(7^(l)+)c/(l2)t + (7,(pw) + /„(7^(2)-))" 

-(-i)"(c/(i2)g_i(r(i)+)[/(i2)t + u_^)(^p(pw) + j^^^(7^(2)-)^ +7'(^^7^(l)-7'i7^7^(l) 



)) 



+ (nlp"^7^(l)-7'i7'(l + r(l)+^)[/(l2)t7^(2)-t(^l + t(l)-t) +h.c.y 



)7^(2)-5_l(r(l)+)7^(2)-t 



where we have defined the auxiliary functions 

' _7^p(pw)^t + 7^7?(<=-) _ 7?i7)7^t a = -l 

^ (-l)'"7^tn*^^7^ + n^vir'' + (-l)"7'i™^7^ a^me {i, 2} 



fain) 



a = -l 



{'y--p(pw)-y-t _ jyipvr) ^ ^p(pw) j-f _|_ p(pw)^j ^ j--p(pw ^ _ 
T^v!r^T - vir^ = f^v!r^f + vi^^^f + f^vi^^^ a = m e {i, 2} 

' jj(12)p(pw)^(12)t _ .p(pw) ^ p(Pw)^(i2)t-^(2)-t^(l)+t 
+7^(l) + 7^(2)-[/(12)p(pw) ^ 7^(l) + 7^(2)-[;(12)p(pw)j^(12)t^(2)-t^(l)+t a ^ -I 

L/(21)t7?,(P-){/(21) _ 7?,(P-) ^ ^^Pw)-^(2)--^(l) + f^(21) 

+C/(2l)t7^(l)+t7^(2)-tp^P-) + c/(2i)t7^(i)+t7^(2)-tp,(pw)7^(2)-7^(i)+^(2i) a = m e {1, 2} 



(85) 
(86) 
(87) 



Equations (83)- (88 1 allows then to explicitly consider two bodies of arbitrary geometries and dielectric properties. 
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in a system characterized by tree possibly different tem- 
peratures Ti, T2 and T3. In order to obtain the expression 



of the force and heat transfer on body 2, in Eq. (85) the 



indexes 1 and 2 must be interchanged, as weU as the in- 
dexes -I- and — . Moreover, in the case of the force the 
overall sign has to be changed. In what follows wc ana- 
lyze such expression for several cases. 



IX. SOME APPLICATIONS 

In this section we are going to perform some applica- 
tions of eqs. (83) and (84). In particular, we are going 



to discuss the force acting on a neutral atom in front of 
a planar slab of finite thickness, as well as the force and 
the heat transfer in the case of two parallel slabs. To 
this aim we will provide the reflection and transmission 
operators associated to an atom and a planar slab. 



VIII. FORCE AND HEAT TRANSFER ON A 
BODY ALONE OUT OF THERMAL 
EQUILIBRIUM 



Before discussing some numerical applications of eqs. 



(83) and (84) for particular choices of bodies 1 and 2, 



we will start by applying our formalism to a simpler con- 
flguration, providing an interesting example of the role 
played by the reflection and transmission operators in 
the calculation of the force and the heat transfer. We 
are now going to consider the problem of a body (called 
body 1) at temperature Ti placed in absence of body 2 in 
the same environment as before having temperature T3. 
In order to obtain the force and the heat transfer in this 
case we can exploit the result (85 1 and impose 72. "^^'^ 
t'(2)± = as well as T2 = T3. We remark that the equi- 



0, 



librium contribution to the force in ( 83 1 goes to zero in 
this limit. After straightforward manipulations, the re- 
sult reads 



.2 — m 



«3 



-L ' rn 



A™(Ti,r3) = (-l)'"+ifiTr{c 
X (-l)"7^(l)+7?^P|")7^(l)+t - 7^(i)-7?(P")7^(i)-t 
+ (_l)™|^7-(i)+-p(p^")7-(i)+t _^ j-(i)+-p(P|"' ^'piP^'^fW+f 

_ ^7-(i)-p(p^™)7-(i)-t _^ f-(i)-p(pw) _^_'p(p:^)f{i)-t^' y 

(89) 



Clearly, this expression is in general different from zero. 
Considering the particular case of thermal equilibrium 
Ti = we see that goes to zero: as expected, at 
thermal equilibrium no force is acting on a body alone in- 
dependently from its geometrical properties, and it does 
not exchange any heat with the environment. On the 
contrary, if Ti T^, the force and the heat transfer are 
linked to the different behavior of reflection and trans- 
mission on the two sides of the body. In particular, if 
the body is symmetric with respect to a plane z = zq, 
it is easy to show that the matrix elements of reflection 
and transmission operators on the two sides cancel each 
other in the case of the force (m =1). This is expected 
for evident reasons of symmetry. Nevertheless, even un- 
der this specific assumption, the heat transfer (m = 2 in 
eq. (89)) still remains different from zero. 



A. Force between an atom and a slab 

Let us start with the case of a neutral atom (body 
2) in front of a slab (body 1) having finite thickness Si. 
The atom has position — (r^,2:^) = (0, 0, z^) (we 
have chosen = in virtue of the cylindrical symme- 
try of the problem with respect to the axis z = 0) with 
ZA > 0, whereas the slab is defined by the two interfaces 
z ~ and z ~ — as shown in figure |4] This configura- 



B 




Za 



z=0 



Figure 4: Geometry of the atom -slab configuration. 

tion is interesting since it implies the presence of a body 
(the atom) not characterized by translational invariance 
and then for which the plane-wave basis is not a natural 
choice. Nevertheless, we will show that the knowledge 
of the atomic scattering operator in this basis, chosen in 
our calculation for convenience, allows us to reproduce 
the known results in some particular limiting cases and 
to give the general expression in presence of three differ- 
ent temperatures Ti, T2 and T3. 

We now first discuss the reflection and transmission 
operators for the slab T?,'^^-'"'" and T^^'''^ ■ For homogeneous 
flat slabs, these operators are diagonal and given by 



(p,k|7^(l)+|p',k') = (27r)2^(k-k')VPip(k,^) 
(p,k|r(i)+|p',k') = (27r)2<5(k-k')V^ip(k,^) 



(90) 



and deflned in terms of the Fresnel reflection and trans- 
mission coefficients modified by the finite thickness 5i 

2 _ p2ifc2i<5i 



1 - (k, u;)e2»fe-i'5i 



Tip(k,w) = 



tip{k,uj)tip{k,uj)e 



(91) 



(k,a;)e 



2ikziSi 
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In these definitions we have introduced the z component 
of the K vector inside medium 1 



£1(0;)- 



k2, 



(92) 



the ordinary vacuum-medium Fresnel reflection coeffi- 
cients 



?'l.TE 



^21 



'"l.TM 



£i(a;)fc^ + fc^i 



(93) 



as well as both the vacuum-medium (noted with i) and 
medium-vacuum (noted with t) transmission coefficients 



Cl.TE 



t 



l.TE 



kz + kz\ 
2fc2i 



'1,TM = 



2^£l(tj)fc^ 



t 



1,TM 



_ 2^ ex(ijj)kz i 



i{uj)kz + kzi ' 



The other ingredient of our calculation is represented 
by the atomic scattering operators. As discussed in [i5] . 
these operators can be deduced, in dipole approximation, 
starting from the description of the atom as an induced 
dipole d{uj) = a(a;)E(R^, w) proportional to the com- 
ponent of the electric field at frequency uj calculated at 
the atomic position Ha- The proportionality factor coin- 
cides with the atomic dynamical polarizability a{uj). The 
field radiated by the induced dipole can thus be written 
analytically as a function of any incoming field: this pro- 
duced field has then to be decomposed in plane waves. 
The expression of the outgoing amplitudes as a function 
of the incoming ones provides the explicit expression of 
the atomic reflection and transmission operators. They 
read (for = -f , — ) 



(k,p|7el(.)|k',p')-^^^^^ 



2eoc'^k, 

X exp[i(k' — k) • r^] exp[—i(j)(kz + k'^)zj<^ 



(e^(k,^).e/(k',..)) 
-i(j>(kz + k'JzA] 
e^(k,c.).6^,(k',c.) 



{k,p\TX{u;)\k',p'] 

2eoc K 

X exp[i(k' — k) • r^] exp[~i(f>{kz — k'^)zA]- 



(95) 



We remark that the we provided the modified atomic 
transmission operator T4 (as a matter of fact, it clearly 
goes to zero in absence of the atom) and that both op- 
erators are not diagonal with respect to the wavevector 
k and the polarization p, as a result of the lack of trans- 
lational invariance on the x — y plane. Moreover, since 
we have attributed a temperature T2 to the atom, the 
atomic polarizability a{uj) must be the one associated to 
a thermal state at the same temperature. 

We are now ready to calculate the equilibrium and 
non-equilibrium force on the atom. Coherently with the 
dipole approximation, we have to keep only the leading- 
order terms in these expressions with respect to the 
atomic polarizability, and thus to its scattering opera- 
tors (95). Moreover, the appropriate changes have to be 



made in eq. ( 85 ) , considering that we are in this case cal- 
culating the force on the body 2. As shown in [45], this 
procedure leads to the expression of the force on the atom 
at thermal equilibrium deduced using several different in- 
dependent approaches. Focusing on the nonequilibrium 
contribution, after some simple algebraic manipulations, 
we obtain 



A2(Ti,T2,r3 



Im 



duj uj'^ a{uj) 



ni3 / " dkk{\pip\^ + \Tip\^ -1) 



+ / dkk{e+ ■ e;) (nsipip e^^'^^'^ + n^zpl^e-^'''''^) 
Jo 



(96) 

where the dependence on the variables uj and k of all the 
quantities inside the integral is kept implicit. The first 
term in the square bracket in eq. (|96| does not depend 
on the atomic position za and it was already identified in 
|9j. On the contrary, the second and the third terms do 
depend on the atom-slab distance za, but they come from 
different regions of the spectrum: the former results from 
propagative waves only, the latter from the evanescent 
sector. As a check of coherence with previous results, 
we have reobtained the expression deduced in [3] using 
a different approach. To this aim, we assumed that the 
atom occupies its ground state (T2 = OK) and that the 
slab and environmental temperatures Ti and T3 are such 
that no atomic excitation is possible: this corresponds to 
the replacement of the frequency-dependent dynamical 
polarizability a{uj) with its static value a{Q). 



B. Force between two slabs 

The case of two parallel homogeneous dielectric slabs 
of finite thickness will be now examined. This configura- 
tion, already studied in [lOj in the case of infinite thick- 
ness, shows the advantage of keeping the translational 
symmetry, making all the scattering operators diagonal 
in the (k,p) basis, allowing at the same time to study 
the effect of the environmental temperature, in virtue of 
the finite thickness of the slabs. Let us assume that the 
slab i [i = \,2) has thickness Si and call d the distance 
between the two slabs: in particular, the slab 1 occupies 
the region —5i < z < (as in the atom-slab configura- 
tion described in sec. IX A^ whereas the slab 2 coincides 
with d < z < d + §2, as shown in figure [5] The reflec- 
tion and transmission operators 7?.^^^+ and 7"^^^^ for the 
slab 1 are still given by eq. (90 1. The operators TZ^^^~ 
and T'-'^-'" associated to the 



eft side of body 1 as well 
as all the scattering operators of body 2 can be obtained 
by solving the problem of the behavior of the scattering 
operators with respect to changes of frame of reference. 
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i.e. the force per unit of surface, is a finite quantity. Nev- 
ertheless, the simple analysis of the symmetrized average 
of the zz component of the stress tensor ([T7| shows us 



Figure 5: Geometry of the slab-slab configuration. 



that this quantity is in this case finite and independent 
on r, coherently with the translational invariance. More- 
over, the result for the average value of T^^, which means 
the force per unit of area, is the same we would get by 
using the formula ( [74| for the force and neglecting the 
divergent term (27r)^(0). 

We are thus now ready to give the explicit expression 
of the pressure acting on slab 1 given by eq. (74 1 after 
neglecting the divergent term. As for the equilibrium 
contribution at temperature T, it is given by 



discussed in appendix The result is that the matrix 
elements of T*-^-*^ coincide with the ones of 7"*-^-'^ given 
by eq. (901, while the interchange of 1 and 2 provides 
directly the elements of the transmission operator T'-'^'". 
As far as the reflection operators are concerned we have 

(p,k|7^(l)-b',k') = (2^)2<5(k-k')VPip(k,^)e-'^'^'^ 

(p,k|7^(2)-|p^k') = (27r)2<5(k-k')VP2p(k,w)e2^'=^'^. 

(97) 

Before moving to the explicit calculation of the force 
we remark that in this geometrical configuration the 
matrix clement of any scattering operator between the 
states |k,p) and |k',p') is proportional to the Dirac delta 
(27r)2j(k - k'), as evident from eqs. ^ and This 



property reflects indeed the translational symmetry with 
respect to x and y characterizing this system. As a con- 
sequence, the total force (74) acting on body 1 is pro- 
portional to (27r)^(5(0) and then formally divergent. This 
happens since we are calculating the total force on slab 
1, which is by definition infinite, while the force density. 



X iV(a;,r) 



(27r)2 ^ 



(98) 



where 



Dp = 1 



„2ik,d 



PipP2pe""-"' (99) 
and the quantities pip, p2p and Dp implicitly depend on 
uj and k. We now turn to the the non-equilibrium con- 
tribution AFi^(Ti, Tz, T3) appearing in eq. ([74|. This 
gives in this case the non-equilibrium pressure 



APi,(Ti,r2,T3) = A(°")(Ti) - A(™)(r2) 



bJp")(T2) 



B 



(pw) 



(T3)-i3^P")(T, 



where we have defined 



B. 



(pw) 



(r3)-B^p*)(ri) 



(100) 



b['""\t) 



B^r\T) 



b^I'"\t) 



h 

2^ 



,Im(pipPSp)^_2rfImfc. 



\Dp? 



h 
h 
h 



.+00 .+00 J 
y duj I dk klmkz n{uj,T)- 

p -'0 h 



duj / dkkkzn{uj,T)- 



\Dp? 



+00 



y I doj dkkkzn{uj,T) 
p Jo Jq 

y dui dkkkzn{uj,T) 



\rip\'{l + \p2p\Hl~\rip\')) 



g^(l-KKp-|ri,n-l 



\ Up 



(101) 



We start noticing that the last term of the last line rium force (100) a contribution 2(j(T^ — T2)/3c where 



can be explicitly integrated. It gives to the noncquilib- a = tt^ kg / 6Qc^ fV^ : this term is the well-known Stefan 
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Boltzmann radiation pressure. We have verified that in 
the hmit of infinite thickness, corresponding to rip, T2p — > 
and the replacement of pip and p2p with the ordinary 
Fresnel coefficients, we analytically reobtain the results 
already deduced in [TUl [H] . 

We have then numerically evaluated, using eqs. (98) 
and (100) arranged as in eq. (83), the total pressure act- 



ing on a 2/im thick slab 1, made of fused silica, in front of 
a 1000/zm thick slab 2, made of silicon. The optical data 
for the two materials are taken from [38]. We have con- 
sidered different sets of temperatures (ri,T2,T3). The 
results are shown in figure [6] and [7j 
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Figure 6: (color online) Slab-slab geometry (see Sec 
Pressure (100 1 acting on a iSi = 2/im thick slab (body 1 
fused silica) parallel to a ^2 = 1000 /im thick slab (body 2, 
silicon). Lines; equilibrium pressures at T = K (black solid), 
300 K (blue dashed), 600 K (red dash-dotted). Symbols: non- 
equilibrium pressures, Ta = OK (blue circles), 300 K (green 
diamonds) , 600 K (magenta plus) , with Ti = 300 K and T2 = 
OK in (a), Ti = OK and T2 = 300 K in (b), Ti ^ T2 = 300 K 
in (c). 



In figure [6] a wide range of distances, from 1 to 
10 jjLTn, has been considered, for different equilibrium and 
nonequilibrium thermal configurations. In particular, 
(a), (b) and (c) correspond to three different choices 
of the slab temperatures Ti and T2 (see caption of fig- 
ure |6] for details). For each case we have represented 
the nonequilibrium pressure corresponding to the values 
of the environmental temperature T-^ = 0,300, 600 K, as 
well as the equilibrium pressure at the same three tem- 
peratures. We note that the transition from an equilib- 
rium to a nonequilibrium configuration can dramatically 
change both the qualitative and quantitative behavior of 
the interaction. It is worth noting that, even for fixed 
values of Ti , T2 and the slab-slab distance d, the value 
of T3 may significantly affect the value of the pressure. 
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Figure 7: (color online) Zoom of figure|6]in linear scales, with 
the same conventions. Here the change of sign of the force 
clearly appears. 



even by orders of magnitude. As a consequence, the envi- 
ronmental temperature can be remarkably considered as 
an efficient tool to tune the interaction. This feature is 
equally present in the case described in figure [6]jc), where 
Ti equals T2. This underlines that even in experiments 
devoted to the measure of the force at thermal equilib- 
rium the environmental temperature should be carefully 
controlled. 

All these effects prove to be even more spectacular by 
looking at figure [7] where linear scales are employed. In- 
deed, in the case of T3 = OK the pressure becomes ex- 
actly zero at a given distance around 6 /im, and repulsive 
for larger distances. The appearance of repulsive inter- 
actions with nonequilibrium systems has been previously 
showed only for microscopic bodies, and in particular for 
the atom-surface interaction [TT] . Moreover, the possibil- 
ity of drastically reducing the Casimir-Lifshitz force may 
be useful in investigations of hypothetical smaller forces 
of different origins (39n43j . 



C. Heat transfer between two slabs 

In analogy with the force, we have chosen the same 
slab-slab configuration to provide a numerical applica- 
tion of the eq. (84) giving the heat transfer on body 1 
for any choice of Ti, T2 and T3. 



This case was already 
studied in [25j , where the infiuence of the environmental 
temperature T3 was not considered. Also in this case, 
the result is expressed per unit of surface: we then ob- 
tain the energy hi absorbed per unit of surface and per 
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unit of time by the slab 1. Its analytic expression, using 
the same formalism of sec. |IXB[ reads 



where we have defined 



b: 



(pw) 



B^r\T2) 



(r3)-6^P")(Ti) 
W 



B^r\T) 



B^r\T) 



27r2 



47r2 



47r2 



E 
E 
E 



-00 Z' + 'OO 



47r2 



E 



duj / dkkioniuj^T)- 



Re(pipP2p) - R.e(pipp;p] 



"'0 

+00 







+00 







duj / dk k u! n{(jj , T) 



|ri,P(l-|p2pP(l-|Ti,n) 



(103) 



* 2 



/ dkkioniLO,Ty^il-\p,p\'-\np\^). 



We have numerically evaluated the heat transfer ( |102p 
on body 1 for different set of temperatures (Ti, T2, Ta). 
In particular, the cases (a), (b) and (c) of figure [s] cor- 
respond to three different choices of the slab temper- 
atures Ti and T2 (see caption for details). For each 
case we have represented the heat transfer correspond- 
ing to the values of the environmental temperature T3 = 
0, 300, 400, 500, 600 K. As in the case of pressure, the role 
of the environmental radiation for heat transfer is partic- 
ularly interesting. In figure [8] the heat transfer hi shows 
an oscillating behavior with an amplitude increasing with 
the temperature T3. As far as the positions of minima 
and maxima are concerned, they are almost insensitive 
to the slab thicknesses 5i and S2 and to the three temper- 
atures Ti , T2 and T3 , being on the contrary connected to 
the dielectric properties of the two bodies. Furthermore, 
these oscillations originate from the propagative sector, 
as evident from the analysis of case (c): as a matter 
of fact, in this configuration, where Ti and T2 coincide, 
eq. ( 102 ) contains only contributions of pure propaga- 



tive nature {B2 and B^). These oscillations were already 
theoretically studied in [T] . Another interesting property 
emerging from figure |8|^a) is the occurrence of a change 
of sign in the heat transfer. Focusing on the blue crosses, 
corresponding to {Ti,T2,T3) = (300 K, OK, 400 K), we 
observe that at large separations, where the propagative 
waves play a dominant role, the heat transfer is posi- 
tive, i.e. the slab 1 absorbs energy. At smaller sepa- 
rations, and in particular for distances of the order of 
1 pm the heat transfer changes sign, becoming negative 
(i.e. the slab 1 radiates energy). This can be understood 



in terms of the evanescent-wave coupling between body 
1 and body 2, which is at zero temperature. We also 
note that the higher is the value of T3, the smaller is the 
distance at which the change of sign occurs. From the 
figure, we deduce that for T3 = 500, 600 K this transition 
happens at distances below 0.5 fiin. 



X. CONCLUSIONS 

We presented a systematic derivation of the radiative 
heat transfer and the Casimir-Lifshitz force between two 
arbitrary bodies. We have first expressed the correlating 
functions of the electromagnetic field in any region of the 
system as a function of the scattering operators of each 
body. This result has been used to provide to closed- 
form unified analytic expression of the heat transfer and 
the force. This expression fully takes into account the 
interaction between bodies of finite size, any shape and 
any temperature, as well as the presence of a thermal 
external radiation coming from the environment. 

We applied this theory to two simple but instructive 
examples: an atom in front of a slab and a couple of 
parallel slabs. The former configuration is an example 
of non translationally-invariant system, and generalizes 
previous results out of thermal equilibrium. The lat- 
ter represents the simplest geometrical configuration in 
which the effects of finite size and external temperature 
can be quantitatively analyzed. As far as the force is 
concerned, we observed that the environmental temper- 
ature can substantially tune the interaction, eventually 



20 



0.2 
0.1 


^ -0.1 

CM 

B 

« 0.2 
0.15 

0-1 
0.05 

X 



(a) . 

• 












♦ ♦♦♦♦ 




" + +* 

;=°°° 


+ + + + + + + + + *■' 


+ + + + + 


□□□□□□□□□□□□□□□□□□□□□□□□□□□□1 













° + * 




° "■■■■□=,=□„„ 


+ + + + + + + + + ++ + + ++ + 4-+ + + + + + + + ++ + + ++ + ++ + - 




4 6 
separation (^jm) 



Figure 8: (color online) Slab-slab geometry (see Sec. IX C I 
Radiative heat transfer per unit of surface ( 102 1 on a Si = 



2 fim thick slab (body 1, fused silica) parallel to a ^2 = 
1000 /im thick slab (body 2, silicon). The temperatures of 
the slabs are Ti = 300 K and r2 = OK in (a), Ti = OK 
and T2 = 300 K in (b), Ti = T2 = BOOK in (c). Symbols: 
Tg = OK (black squares), 300 K (green plus), 400 K (blue 
crosses), 500 K (brown diamonds), 600 K (red circles). 



producing a repulsive force. As for the heat transfer, it 
shows an oscillatory behavior with respect to distance 
connected to the dielectric properties of the two slabs 
and whose amplitude increases with the temperature. 
Moreover, we found that some given choices of the en- 
vironmental temperature are able to produce a heat flux 
whose sign changes as a function of distance. 

This study shows the interest of nonequilibrium con- 
figurations, mainly consisting in a strong tunability of 
force and heat transfer. It would be thus interesting to 
apply our results to other geometrical configurations and 
to test them experimentally. 
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Appendix A: Transformations of scattering 
operators under translations 

We will discuss here the problem of the transformation 
of the scattering (both reflection and transmission) op- 
erators with respect to translations. Let us suppose to 
have a first frame of reference O and a second one O ob- 



tained from the first one by a translation of Rs- For the 
purpose of this section it is useful to gather the refiec- 
tion and transmission operators TZ^ and in a single 
operator iS*^*^ connecting the outgoing modes propagat- 
ing in direction (j) to the incoming modes propagating in 
direction cj)' . As a consequence, we have the following 
identification 



S++ = 

S-+ = n- 



(Al) 



s+- = n 
s— = r 

From the definition of the electric field in O 
E(R,a;)=^ J ^ exp(zK^ • R)e^(k, a.)ii;,^(k, ..) 



<t>,p 

E 

0,P 



exp[zK^.(R-Rs)]enk,a;) 



exp(zK"^-R5)£;;J(k,( 



(A2) 



we deduce that the amplitude i?^(k,w) in the new frame 
of reference O equals 



(A3) 



E^ik,Lu) = exp(iK* • Rs)i?p*(k,w). 
From this properties we deduce 

^^°"*)*(k,cj) = exp(iK* • Rs)4°"')^(k,cj) 

= exp(*K^.R5)5] / ^{p,k\S^^'\p',k') 
p' 

x4")^'(k',..) 

= exp(*K^.R5)5] / ^(p,k|5^^'b',k') 
xexp(-*K>'.R5)^^!"^^'(k',c.) 



and finally the link between the matrix element of the 
scattering operators in the two frames of reference 



(A4) 



(p, k\S^^' \p\ k') = exp[i(K* - k'*') • Rs)] 
X (p,k|5*^'b',k'). 



(A5) 



In the particular and important case of translation along 
the z axis, assuming that the origin of O has coordinates 
(0, 0, d) with respect to O, we have 

(p,k|^+|p',k') = eMKkz + fc;)d](p,k|7^+b',k') 

(p, k|7^- k') = exp[-i{k, + K)d] {p, k\n- \p', k') 
(p,k|r+b',k') = eM^ikz - K)d\{p,^r+\p'y) 

{p,k\f-\p',k') = exp[-z(fc, - K)d]{pMT-\p',k'). 

(A6) 

Apart from their general theoretical interest, these rela- 



tions are used in Sec. IX in order to deduce any reflection 
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and transmission operator associated to a planar slab as 
a function of the ordinary Fresnel coefficient (modified to 
take into account the finite thickness), usually calculated 
assuming that the interface coincides with the surface 
z = 0. 



Appendix B: Reciprocity relations of scattering 
operators 

The matrix elements of the scattering operators are 
not all mutually independent. By exploiting some prop- 
erties of the electromagnetic field, it is in fact possible to 
deduce some relations connecting these elements. This 
is the case, for example, of the reciprocity relations pre- 
sented in |44j . In this appendix, we derive and express 
these relations using our field decomposition and nota- 
tion. In [33], the authors start their derivation of the 
reciprocity relations from Lorentz's reciprocity theorem 
in presence of sources. To formulate this theorem, we 
start by supposing to have, in presence of a given body, 
a dipole pi (P2) in position Ri (R2) and oscillating at 
frequency w. Each of these dipoles produces an electro- 
magnetic field which is then scattered by the body and 
reaches the other dipole. Lorentz's reciprocity theorem 
can then formulated by imposing that 



Pi-E2(Ri) = P2-Ei(R2) 



(Bl) 



where Ei(R2) (E2(Ri)) is the result of the scattering 
on the body of the field produced by the dipole pi (P2), 
then calculated at the position of dipole P2 (pi)- 

In order to fix the notation, let us suppose that the 
dipole pi is on the side (j) of the body, while p2 is on the 
side 4>' . In order to calculate the field Ei(R2), we need to 
take the component of the field emitted by pi propagat- 
ing in direction —cj), then apply the operator S'^ and 
calculate the resulting field in position R2. The compo- 
nent at frequency w of the field produced by a dipole p 
in position Rp reads 



E(R,a;) 



-Vr X Vr X 



Rd 



(B2) 



where Rd = |Rd| = |R — Rp| and Vr represents the 
gradient with respect to R. As discussed in [45] , the pas- 
sage from the spherical wave emitted by the dipole and 
our angular spectrum representation can be performed 
by using the Weyl representation ^46]: the resulting field 
propagating in the direction from position Rp reads, 
after straightforward algebraic manipulations, 



E'^(R,a;) 



lU! 

2^ 



(e^(k,L.).p) exp[iK^. (R-Rp)]. 



(B3) 



We can then simply deduce the expression 



Ei(R2) 



E 



pp 



X 



^(pi.6/(k',a.))(k,p|5^'-^|k',p') 
xexp[i(K'^' •R2-K'-*'-Ri)]. 
and analogously 

■2 r r d2k' 



(B4) 



E2(Ri) = 



E 



pp 



2eoc2 J (27r)2 J (27r)2 



e(-k',c.) 



X l(p2 • e;^'(-k,c.))(-k',p'|5^-*'| -k,p) 
xexp[i(K'*' •R2-K'-^-Ri)]. 
where we have performed the change of variables 
(k,k',_p,p') — > (-k',-k,p',p). 



(B5) 



(B6) 



By imposing the condition (Bl I we get 



K (pi • V^k', io)) (p2 • < (k, uj)) (k, p\S^' 
= fc;(pi.e^,(-k',c.)) (p2-ep^'(-k,a.) 
X (-k',p'|5^-^'|-k,p). 



(B7) 



from which we deduce, by using the properties of the 
polarization unit vectors ([8|, the final relation 

fc,(k,p|5*''-^|k',p') - fcU-ir+P'(-k',p'|5*-^'| - k,p). 

(B8) 

By choosing the four possible values of the couple (0, </)') 
we obtain the relations 

fc,(k,p|r±|k',y) = fcu-if +p'(-k',p'|r^l - k,p) 
fc,(k,p|7^±|k',p') = fc;(-if +f'(-k^p'|7^±| - k,p). 

(B9) 

These equations clearly show that the matrix elements of 
each transmission operator are connected to elements of 
the specular operator, while the matrix elements of each 
reflection operator are not independent. The relations 
(|B9|) contribute to the derivation presented in appendix 



Appendix C: Green function and scattering 
operators 

In this appendix we are going to derive the relation be- 
tween the Green function in the presence of a single body 
and its reflection and transmission operators. This fea- 
ture is a main point in the calculation given in appendix 
Id] In order to derive this connection we will start from 
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the definition of the Green function. Suppose to have a 
dipole electric moment p located at R' oscillating at fre- 
quency w and thus producing an electric field oscillating 
at the same frequency proportional to the components of 
the dipole moment itself. The component Gij(R, R') at 
frequency lu of the Green function can be interpreted as 
the part of component i of the total electric field at the 
point R, namely (R) , proportional to the component 
J of the dipole moment, divided by pj. Of course, in our 
description of the electromagnetic field in the presence 
of scatterers, the field directly produced by the dipole 
p will result in reflection and transmission: as a conse- 
quence, the Green function will prove to be linked to the 
scattering operators H,"^ and T"^ ■ 



It is important at this point to remind that our choice 
of mode decomposition of the field naturally introduces 
a left and a right side for a given body. Thus, we will 
separately discuss the cases in which the arguments R 
and R' appearing in the Green function are either on 
the same side or on opposite sides of the body. Let us 
suppose first that the two points R and R' are located 
on the same side of the body. In this case, the field 
directly produced by the dipole at the point R' will be 
directly observed in R. Moreover, this field will produce 
a reflected field defined in the same region (/>, and conse- 
quently observed in R as well. We argue then that for R 
and R' on the same side of the body the Green function 
can be expressed as a sum of two terms, a free one inde- 
pendent on the scattering operators, and a reflected one 
proportional to TZ"^ . Analogously, if the flrst argument 
R of the Green function is located on the side </>, while 
R' is located in the —(p region, the Green function will 
be made up of a unique transmitted term, proportional 
toT"^. 



R and R' on the side (/> of the body reads 

G,, (R, R', c^) = G|f (R, R, c^) + G,[f (R, H' ,uj), 



E 



d^k 



exp[ik • (r — r')] 



X l [e{z - z') (e+ (k, uj))^ (e+ (k, w)) ^ exp [ik, [z - z')] 
+ 9{z' - z) {e^ (k, w)^ (k, w)^ exp [ik^ {z' - z)] 



Gif(R,R',a.) 



E 



d^k 



d^k' 



2eoc2 ^7 (27r)2 J (27r)2 
exp[i(k . r - k' . r')] ^ (e^(k, a.)) ^ (e;/(k', c.) 
e-xp[i4>{kzZ + k'^z)\ (k,p|7^'^|k',p'). 



(CI) 



On the contrary, if R (R') is located on the side (f) {—<j>) 
of the body, we observe a fleld, propagating in the 
direction, resulting from the transmission described by 
the operator T"^ of the component of the fleld emitted 
by the dipole propagating in direction (j) as well. As a 
consequence, in this case we conclude 



G,,(R,R',a.) - ^ , ^^^^2 



E/|;;;^(s^(^-)).^'(^» 

pf 

ex.p[i<p{k;,z - k'^z')] (k,p|T'^|k',p'). 



(G2) 



The expressions (CI) and (C2) describe, for any position 



of the points R and R' with respect to the body, the 
connection between the Green function and the scatter- 
ing operators. 



Appendix D: Correlators of the field emitted by the 
bodies 



In order to make this description analytic we need the 
explicit expression of the dipole field propagating in di- 
rection (j> given in appendix [Bj As far as the free contri- 
bution is concerned (existing if R and R' are located on 
the same side of the body), if R is on the right (left) side 
of R', the Green function will contain the component of 
the field emitted by the dipole propagating toward the 
right (left). Let us now suppose that both R and R' are 
on the same side (j) of the body. In this case, apart from 
the direct contribution we have just discussed, the field 
contains the component propagating in direction 0, re- 
sulting from the reflection by means of the operator TZ'^ 
of the dipole field propagating in the opposite direction 
—(j). Then, the ij component of the Green function for 



As discussed in sec. |VH when a body is at ther- 
mal equilibrium with the environment at temperature 
T, the fluctuation-dissipation theorem gives complete 
knowledge of the correlation functions of the total field, 
resulting from the environmental field, the one emitted 
by the body, and the result of scattering processes. Nev- 
ertheless, even at thermal equilibrium the correlation 
functions describing the field emitted by the body are 
not straightforward. They were given in eq. (451 for 



two components of the field propagating in the same di- 
rection, whilst ( 46 1 gives the corresponding quantity for 



counterpropagating components. The derivation of eqs. 
( 45 1 and ( 46 ) is the main scope of this appendix. 



Let us consider a body at thermal equilibrium at tem- 
perature T and described by the scattering operator TZ^ 
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and ■ The first step of our calculation is to write the 
expression of the total field both on the left and on the 
right side of the body. Their components at frequency uj 
are given by 



+ E(''<=)'*(R,a;) + E("')'*(R,a;). 



(Dl) 



In this expression the superscript for the total field in 
the l.h.s. refers to the region we are considering, while 
for the four fields in the r.h.s. it corresponds, as usual 
in the rest of the paper, to the direction of propagation. 
Eq. (Dl I tells us that the total field in region cf) contains 



first the environmental field propagating in direction —0 
(i.e. toward the body) and of course the field eC^)"*^ emit- 
ted by the body itself and propagating in direction ip. 
Moreover, the environmental field E^'^'™'"'^ produces a 
reflected field E^"'")'^ connected to E^™^'"'* by the opera- 
tor TZ'^. Finally, in region </> there is a field E'^'''^'^ resulting 
from the transmission of the environmental field E^'^"^)''^ 
(existing in the other region —<j>) by the operator T"^. 
We are looking for correlation functions such as 



(4^)nk,^)i?^'^*'^(k',a.')) 



sym 



(D2) 



of the field emitted by the body, both for <j> — (p' and for 
(j) ^ (j)' . Let us start by the case of copropagating com- 



ponents of the field, namely by the case (f> = (f>' . It is first 
useful to calculate the correlation function of the i and j 
components of eq. (Dl), for two couples of coordinates 



(R, u) and (R', u'). As for the l.h.s., the result is directly 
given by the fluctuation-dissipation theorem stated in eq. 
(39). The discussion of the r.h.s. requires a more accu- 
rate analysis. We first point out that we already know 
the correlation function characterizing the environmental 
field: it is given by eq. ( 43 ) and it also tells us that two 



counterpropagating components of this field are uncor- 
related. As a consequence, E*^""^)""^ is correlated (apart 
from with itself) with the reflected field E^''')'^, but not 
with the transmitted field E^*''^"^, coming from the coun- 
terpropagating component of the environmental field. As 
a consequence E''"')"^ is correlated only with itself: the 
same holds for the field E*^^)*^, since the field produced 
by the body alone is of course uncorrelated with every- 
thing which results from the environment. We now ob- 
serve that, since the reflected and transmitted field are 
connected by means of scattering operators to the envi- 
ronmental field, whose correlation functions are known, 
the only unknown remaining is the correlation function 



(D3) 



strictly connected to the one we are looking for. In par- 
ticular we obtain the expression 



L = R, 

R ^ (£;f"")-*(R,cj)i;]^"")-^^(R',c^')) + (i;f^)^(R,c^)£;j'-'=)'^^(R',cj')) + (^|™'^"'^(R,c^)£^f'^^^(R',t^')> 

+ (£;f'=)^(R,w)£;f""'-'^^(R',c^')) + {El'''>'^{R,u)Ef''>'^\-R',Lj')) + ^"^(R, c^)£;f ^^^(R', c^')> 



where the subscript (sym) has been dropped for simplicity. We will start by working on the r.h.s. term R. In this 
case, expression all the fields by their Fourier decomposition ^ (all for a given value of (j)), it is easy to see that all 
the correlation functions will contain a double integral on k and k' as well as the factor exp[z(k • r — k' • r')], which 
we will drop here, in virtue of the fact that it will also appear in the l.h.s. L. Moreover we will also drop the factor 
27r(5(a; — uj'), appearing in both L and R. With these conventions and using the definition of the scattering operators 
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given in Sec. [V]we have 





)) = ^ exp[— i(/)(fczZ — 
pp' 


k':z)](e-'^{k,uj] 


) fe-'^(k',t^)V(k,p|C(3)|k',p'), 




)) = eyiY>[i4>{kzZ — k 
pp' 




^et, (k', w)) * {\i,p\n^c^^'^n^^ |k', p' 




)) = exp[— i(/)(fczZ 4- 




) fet(k^c^)y(k,p|C(3)7^'^t|k',p' 




)) — ^ exp[i0(/czZ + k 
pp' 


l*z')](e^(k,a;)) 


V/(k',c.))*(k,p|7^*c(3)|k^p'), 


(£;f'')^(R,c^)£;f^'^^(R',L.' 


)) = exp[i(j>{kzZ — k 
pp' 


l*z')](ep(k,w)) . 


'6^,(k',a;))*(k,p|r*C(3)r^t|k>' 


(^f^^R,cj)i;f^'^^(R',cj' 


)) — exp[i0(A:2Z — /; 
pp' 


l*z')](e^(k,a;))^ 


^e^,(k',a.))*(k,p|C^^|k',p'). 



(D5) 



Wc now need to calculate the l.h.s by using the of G-°^ finally takes the form 



IniGg'(R,R',w) 



X 6{uj — ck) 



1 



fluctuation-dissipation theorem (|39| and inserting the ex- 
plicit expression found in appendix [C] of the Green func- 
tion as a function of the scattering operators. Since in 
this case we are calculating the correlator of E'f°*^^'^ with 
^jtot)0^ we are using the Green function calculated at 
positions R and R' being on the same side (f) of the 
body. We will thus make use of eq. (CI) to connect 
Green function to scattering operators. In particular, the 
fluctuation-dissipation theorem tells us that we need to 
calculate the imaginary part of the Green function (CI I. 
Starting from the free term cf^, and inserting a term 
(27r)^(5(k — k') integrated over k' as well as a 5ppi in or- 
der to have the desired factor exp [i(k ■ r — k' • r')] and 

the same structure of the terms in (D5| (we remark here Observing from eq. (431 that 
that the factor 2'k5{lj — lj') is already in the fluctuation- 
dissipation theorem (|39P), we obtain 



^V(2^)''^(k-k') 



(D7) 



: (k, u))^ (e+ (k', Lo))^ exp {k,z - k',z')] 
(k,cj)j _(^ep,(k',a;)^ exp[i(fc^z' - k^z)] 



ImGi^")(R,R',a.) = — ^ V (2^)''^(k " k') 

pp' 

X {0(z-z')[^(ep(k,c.))^(e+(k',c.))^exp[z(fc,z-fc;z')] 



(k,p|G(3)|k',p') = -^7V(c.,r)V(2^)''5(k-k') 



X 9{lo — ck)- 



(D8) 



+ 9{z - z 
1 



r 1 



+ ^ [e- (-k, ^)) * (cp, (-k', c.)) * exp [-i{k':z' - kU)]]} 



and taking back the factor -N{uj,T) in cq. (39) we con- 
exp[i(fczZ — kzz)] elude that the term (D7) cancels the first term in (D5) 



and gives in the l.h.s. a contribution 



(D6) 



J2 ewMkzZ - k':z')] (e^(k, lo)) ^ (e^, (k', uj) 
X (k,p|C(3)|k',p')- 



pp' 



(D9) 



where for the terms obtained by complex conjugation we 
made a change of variable from k to — k. Making use 
of the properties of the polarization unit vectors, it is 
easy to show that each term multiplying a 9 step func- 
tion is zero in the evanescent sector of k (and then k'). 

After simple algebraic manipulations the imaginary part We are now left with the calculation of the imaginary 
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part of G^f ^ . We have 



appendix [B] to conclude that 



ImClf (R,R',cj) = 



E 



^ (e^(k,c.))^(e/(k',c.))^ exp[z</.(fc,z + fc^z')] 
X (k,p|7^*|k^p') 

i;^(e;nk,^))*(e^,(k',^))* expH0(fc:z + 

X (-l)f+P'(-k^p'|7^*t|_k,p)". 

(DIO) 



where the second term was obtained by the change of 
variables (k, k') — > (— k, — k') and using the properties 
of the polarization unit vectors ([s]). Starting from the 
first term we note that its part which is propagative in k' 
exactly cancels the fourth term in eq. (D5|. Observing 
that in the evanescent sector we have k* = —k'^ and 
using again ^ we are left from the first term with a 
contribution 



4^ 



PP' - (Dll) 

X exp[t<j>{k,z - k';z')] (k,p|7^^7'i7^|k^p')■ 

As far as the second term is concerned, we exploit the 
reciprocity relations of scattering operators presented in 



1 (_1)P+P'(_k^p'|7^^t|_k,p) = i^(k,p|7^*t|k',p'). 

(D12) 

As a consequence the propagative part with respect to k 
of this second term analogously cancels the third term in 
(D5) and we are left with 



pp' 



(D13) 



X exp[i^{k,z - k';z')] (k,p|7'i7^7^■*^|k^p') 



Having considered all the terms in the equality (D4), we 



remain now with quantities which have all exactly the 
same structure, namely the same sum over p and p' , the 
same z-dependent exponential and the same polarization 
unit vectors. As a consequence we can identify the ma- 
trices whose elements are calculated between (k,p) and 
(k',p') and obtain immediately the equality (45). 



In the case of the correlation function between two 
components of the field emitted by the body propagating 
in two opposite directions the structure of the calculation 
is the same. Now one has to use eq. ( C2 ) instead of (CI) 



in order to connect the Green function to the scattering 
operators. The analytic expression ([46| of the correlation 
function for counterpropagating fields has been already 
given in the paper. 
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